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SUMMARY 


This report describes a ray-trace analysis of aberrations and extreme 
rays of a Cassegrain-type telescope with a tilted secondary mirror. The work 
was motivated by the need to understand the factors limiting image quality and 
to assist in the design of secondary mirrors for three telescopes with oscil- 
lating secondary mirrors (OSM) used at Ames Research Center for high-altitude 
infrared astronomy. The telescopes are a 31-cm-diameter Dall-Kirkham (ellip- 
tical primary, spherical secondary) flown aboard a Lear jet, a 71-cm balloon- 
borne Dall-Kirkham flown on the "AIROscope" gondola, and a 91-cm true 
Cassegrain (parabolic primary, hyperbolic secondary) flown aboard a C-141 jet 
transport. The optics for these telescopes were not designed specifically for 
OSM operation, but all have OSM's and all must be used with various detector 
configurations; therefore, a facility that evaluates the performance of a 
telescope for a given configuration is useful. The analytical expressions are 
summarized and results for the above systems are discussed. Details of the 
calculation and a discussion of the computer program are given in the 
appendices. 


INTRODUCTION 


The requirement of subtracting sky-emission noise in infrared astronomi- 
cal observations has led to the development of the oscillating secondary 
mirror (OSM) technique (ref. 1); the secondary mirror in a Cassegrain configu- 
ration telescope is positioned alternately at one of two positions so that the 
detector alternately views two adjacent areas of the sky. The resulting 
"chopped" detector output is typically processed by a synchronous- amplifier 
(ref. 2) referenced to the OSM motion, thereby providing a dc signal propor- 
tional to the difference in the power received from the two areas of the sky. 
The method has been successful in both airborne and ground-based observations. 
The principal advantages of the technique are (1) the amount of stray radia- 
tion chopped is small and (2) all infrared detection systems to be mounted on 
the telescope are simplified because the chopper is an integral part of the 
telescope. 


♦present address: Physics Dept., Stanford University, Stanford, Calif., 

94305. 



However, some questions arise concerning the quality of the image in OSM 
telescopes. For at least one position of the OSM, the axes of the primary and 
secondary mirrors do not coincide. At wavelengths where thermal emission from 
the telescope is significant (lOOOy i X i 3y), it is desirable to balance the 
power received by the detector for each of the two mirror positions; therefore, 
the secondary mirror is typically aligned so that, for the two positions, its 
axis makes equal and opposite angles with the primary axis. The resulting 
aberrations can restrict the useful field of view of the instrument. Although 
diffraction at infrared wavelengths tends to reduce imaging requirements, good 
spatial resolution of many infrared sources is needed to improve our under- 
standing of them; thus, even in the infrared it is desirable to understand 
aberrations of the optics. For example, these aberrations determine the wave- 
length longward of which diffraction effects become dominant. In addition, 
the quality of the image is an important consideration for multiple-element 
detector arrays used in conjunction with an OSM telescope. 

Also, the possibility of offset guiding by tilting the secondary mirror 
is limited by the resulting aberrations of the image. This application is the 
subject of a calculation by Bottema and Woodruff (ref. 3) , who summarize 
third-order aberration constants for well-known Cassegrain- type telescopes and 
derive the requirements for a configuration with a tilted secondary mirror in 
which the on-axis image of an off-axis object point is free from third-order 
coma. 


Another unusual requirement for an infrared telescope is that the primary 
mirror should be "underfilled," that is, the secondary mirror should be under- 
sized to minimize the amount of radiation from the telescope tube seen by the 
detector. This results in a reduced aperture for which the maximum useful 
value depends on the tilt of the secondary mirror. 

The present report describes a ray-trace analysis of aberrations and 
extreme rays of a Cassegrain- type telescope with a tilted secondary mirror. 
This work was motivated by the need to understand limiting image quality and 
to assist in the design of secondary mirrors for three OSM telescopes used at 
Ames Research Center for high-altitude infrared astronomy. The telescopes are 
a 31-cm-diameter Dall-Kirkham (elliptical primary, spherical secondary) flown 
aboard a Lear jet, a 71-cm balloon-borne Dall-Kirkham flown on the "AIROscope" 
gondola, and a 91-cm true Cassegrain (parabolic primary, hyperbolic secondary) 
flown aboard a C-141 jet transport. The optics for these telescopes were not 
designed specifically for OSM operation, but all have OSM's and all must be 
used with various detector configurations; therefore, a facility that evalu- 
ates the performance of a telescope for a given configuration is useful. 
Results appropriate to these systems are discussed. Details of the calcula- 
tion and discussion of the computer program are given in the appendices. 

The authors gratefully acknowledge the assistance of Mrs. Mai Liu, who 
programmed the plot routine for the line printer, Mr. Gordon Augason for use- 
ful discussions, and Dr. R. T. Jones for the loan of his notes on Cassegrain- 
type telescopes. 
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DESCRIPTION OF CALCULATION 


The calculation which follows obtains the intersection with the image- 
plane of a ray which is incident on the primary at a known radial distance and 
inclination from the primary axis. The secondary surface is assumed to have 
its axis tilted from the primary axis by a known angle 0. The image quality 
is then judged from the distribution of points in the focal plane resulting 
from a uniform distribution of rays having the same inclination to the primary 
axis and filling the primary aperture. 

The geometry and notation are given in figure 1. Ray Li from the 
object strikes the primary at point A, ray L 2 strikes the secondary at 
point B, and ray L 3 intersects the image plane at point D. The secondary 



Figure 1.— Coordinates for the calculation. The origin is at the vertex of the primary mirror, 
whose axis is the z axis. The x axis is parallel to the axis of rotation of the secon- 
dary mirror. 


pivots about point E through an angle 0 which lies in the y-z plane. 
Parameters in the calculation are 

c on- axis distance from center of primary surface to image plane 
Dp diameter of primary mirror 

D s diameter of secondary mirror 

ep eccentricity of primary mirror 

e s eccentricity of secondary mirror 
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Fp rp/2 focal length of primary paraxial sphere 
rp radius of curvature of primary mirror 

r s radius of curvature of secondary mirror 

zl on-axis distance from center of primary surface to image plane 

z s on-axis distance from center of primary surface to untilted secondary 

surface 

y object angle in the y-z plane 

All these quantities except z L are typically positive. Ray Lj makes an 
angle a with the normal to the primary at point A and L 2 makes an angle & 
with the normal to the secondary at point B. 

A few simple formulae are useful. The equivalent f number of the 
system is 


* _ z s ' Z L 

* e 5: 


(i) 


The relationship between the radius of curvature of the secondary (r s ) and the 
focal length of the paraxial primary sphere (Fpj is 


Fp “ z s z s - Z L r s 


( 2 ) 


A typical application of equation (2) is to find the distance to the secondary 
given Fp, z l, and r s ; the solution is 


= 7 {Fp + Z L - r s + f r s 2 + { p p - z l) 2 1 1/2 > 


(3) 


The variation in focal distance for a small change in secondary mirror posi- 
tion is given by 


6 



The underfilled primary diameter is 


D 


pu 


FpDg 


(4) 


(5) 
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and the plate scale (angular displacement of the object in arc-min required 
for a 1-mm displacement of the image) is 


S 


3438’ 

" ^*pu 


( 6 ) 


The angle of tilt of the secondary ( 0 ) determined by the beam angle in the sky 
(y) and the on-axis detector is 


e 


r s F P Y 
r s - c Fp - z s 2 


(7) 


Equations (2) to (7) do not consider the figure of the mirrors, but are good 
approximations in practical cases. 

Appendix A gives a brief analysis of Cassegrain- type telescopes similar 
to that of Jones (ref. 4) with a view to determining the mirror parameters for 
a particular configuration. Appendix B gives some approximate expressions 
useful in estimating effects that tend to overfill the primary. Appendix C 
describes the details of the ray-trace calculation. 

The ray- trace program uses the parameters listed under Description of 
Calculation to compute the secondary tilt angle 9 necessary to center the 
image, the unobscured primary area (underfilled primary area minus secondary 
mirror surface area), the f number of the secondary cone, the plate scale in 
the image plane, and the positions at the primary surface of the two extreme 
rays, which are rays reflected to the edges of the secondary mirror surface in 
the y-z plane. Then the program fills the unobscured aperture with a uni- 
form distribution of parallel rays incident at the specified beam angle y. 
Each ray is labeled with a symbol, of which the first 44 are distinct and the 
remainder are asterisks. The program lists the coordinates of the rays inci- 
dent on the primary and then plots this pattern on the line printer. With 
this incident beam, the program calculates the centered image pattern and 
plots the pattern in the image plane. Alternatively, an arbitrary tilt-angle 
6 can be used, for which the program computes an off-axis image. 

The angular resolution of the telescope is diffraction limited for 
wavelengths greater than some value X c , corresponding to the finite image 
size of a point object. The angular diameter (in arc-min) of the first dark 
ring of the diffraction pattern for wavelength A is 


<P * 


2. 44 (3438 ) A 
Dpu 


( 8 ) 


About 84 percent of the radiation falls within an angular diameter given by 
equation (8) for the diffraction pattern from a circular aperture of diameter 
D d . The program calculates a "diffraction diameter" d of a circular 
aperture that contains about 84 percent of the rays in the image plane. Set- 
ting <j> s Sd for A = A c and using equation (6) , one obtains 
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(9) 


Xc = 2.44 f e 

The program uses equation (9) to estimate X c . 


EXAMPLES 


The parameters for the three telescopes considered are given in table 1. 
All distances are given in millimeters. Values for the eccentricities were 
calculated from the formulae of appendix A. Values of z s> f e , 6, and 0/y 
were determined from equations (3), (1), (4), and (7), respectively. 

Figure 2 shows images calculated for the Lear telescope using the 
parameters from the first row of table 1. The input pattern covers the useful 
area of the primary (fig. 2(a)). The finite size of the centered image from 
an on-axis point object (fig. 2(b)) shows that the telescope is not figured 
precisely for this configuration. That the telescope is well focused is 
indicated by the thorough mixing of the rays from the input pattern. The 
diffraction diameter is 0.06 mm, and 85 percent of the rays lie within a cir- 
cle of this diameter. For this configuration, the telescope is diffraction- 
limited for wavelengths greater than about 3.4p. Figure 2(c) shows the image of 
a point object 5 arc-min off-axis when the secondary is untilted. The image 
appears nearly centered at the origin of the plot, which is displaced in this 
case from the axis of the telescope by 2.81 mm, corresponding to the 5 arc-min 
separation of the source from the telescope axis. Note that the scale size is 
twice that of figure 2(b). Figure 2(d) shows the on-axis image (centered by 
tilting the secondary mirror) of an object 5 arc-min off-axis. The scale is a 
factor of 2 larger than in figure 2(c). The aberrations incurred by centering 
the image of an off-axis object by tilting the secondary are coma- like and 
somewhat worse than those where the secondary is not tilted and the image is 
allowed to fall off-axis. 

Figure 3 shows calculated images for the Lear telescope in a configura- 
tion determined by the parameters of the second example in table 1. The large 
back-focal distance (equal to the primary diameter) of this configuration is 
far from that for which the primary was figured. This accounts for the rather 
large on-axis image of an on-axis point source shown in figure 3(a). Fig- 
ure 3(b) shows the on-axis image of a point source 12 arc-min off-axis. For 
this case, d ~ 0.35 mm, and the telescope is diffraction-limited for wave- 
lengths greater than X c * 23y. The large beam throw in this example does not 
appreciably increase the image size over the on-axis case (fig. 3(a)) because - 
as noted above - the telescope is far from optimized at this large back-focal 
distance. 

Theoretical images for the AIROscope telescope are shown in figure 4. 
Figure 4(a) shows the on-axis case and figure 4(b), the centered image of a 
point source 4 arc-min off-axis. The diffraction diameter of the latter image 
is 0.31 mm, corresponding to = lly. In this example, the two images 
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TABLE 1.- PARAMETERS FOR SAMPLE CALCULATIONS SHOWN IN THE FIGURES. 


a 

D , 

r , 

6 , 

D , 

r . 

e > 


Z T > 

Z c j 

f , 

r 1 

D , 

s. 

6, 

e/v. 

Y j 

Telescope 

V 

p 

p* 

s y 

s’ 

s' 


L 

S 

e 

pu 






mm 

mm 

— 

mm 

mm 

— 

mm 

mm 

min 


mm 

1 /mm 




Lear 

305 

919 

0.82 

78 

325 

0 

20 

-186 

336.1 

6.7 

289 

1.77 

18.7 

2.0 

5.0 

Lear 

305 

919 

.82 

99 

435 

0 

20 

-305 

299.6 

6.0 

288 

1.99 

15.2 

1.5 

12.0 

AIRO scope 

711 

2840 

.89 

124 

615 

0 

18 

-284 

1166 

11.7 

693 

.42 

33.6 

2.9 

4.0 

C-141 

914 

3657 

1.00 

187 

923 

1.35 

10 

-1219 

1435 

14.2 

868 

.28 

47.1 

2.4 

3.0 

C-141 

914 

3657 

1.00 

187 

899 

1.28 

10 

-1752 

1435 

17.0 

868 

.23 

65.4 

2.4 

3.0 


a Used in infrared astronomy at Ames. 
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Figure 2.— Calculations for the Lear jet telescope using the parameters from the first example in table 1. (a) Input pattern, 
(b) Image of an on-axis point object with secondary untilted, (c) Image of a point object 5 arc-min off-axis with secon- 
dary untilted, (d) Same as (c) but with secondary tilted to center the image on axis. 
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differ considerably in size because of the aberrations of the tilted secondary, 
off-axis object configuration considered in figure 4(b). 

Figure 5(a) shows images for the C-141 telescope in the design 
configuration, except that the secondary mirror diameter has been reduced to 
underfill the primary. Note that the design configuration results in good 
image quality for on-axis objects, as expected. The centered image for a 
point source 3 min off-axis is shown in figure 5(b). The diffraction diameter 
is 0.16 mm, which corresponds to A c = 4.6u. 

Figure 6 shows images for the last example in table 1. This back- focal 
distance is longer than the previous example and corresponds to an instrument 
somewhat behind the intended image plane. The secondary in the calculation 
has been refigured for this application, but the primary is as shown in fig- 
ure 5. The quality of the image is comparable to that of the design configura- 
tion (compare figs. 6(a) and 6(b) with figs. 5(a) and 5(b)). Although the 
diffraction diameter is larger for figure 6(b) than for figure 5(b), the 
diffraction- limiting wavelength A c is about the same because of the 
difference in f e for these two systems. 


CONCLUDING REMARKS 


The quality of the image calculated by the program described here is 
somewhat better than the systems represented because of imperfections in the 
optics. Image quality similar to that calculated here has been observed with 
the Lear jet telescope and the C-141 telescope. 

The advantage of the ray tracing analysis over an approximate analytical 
treatment is that ray tracing includes the effects of all aberrations simul- 
taneously and yields a quantitative distribution of the energy in the image. 
This permits a simple comparison with the diffraction distribution and pro- 
vides experimenters with a convenient summary of the anticipated optimum 
telescope performance for a given configuration. 

The most serious type of aberration encountered - third-order coma - can 
be reduced using the optics design of Bottema and Woodruff (ref. 3) and should 
probably be done when one designs a telescope for use with an OSM. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., October 18, 1974 
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Figure 5.— Images for the C-141 telescope at the nominal bent -Cassegrain focal plane, fa) Focused 
on-axis image of an on-axis point source, (b) On-axis image of a point source 3 arc-min off-axis. 
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Figure 6»— Images for the C-141 telescope at the n backplate 1f bent-Cassegrain focal plane for the 
parameters of the last example in table 1. (a) Focused on-axis image of an on-axis point 

source, (b) On-axis image of a point source 3 arc-min off-axis. 
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APPENDIX A 


ANALYTICAL ANALYSIS OF CASSEGRAIN-TYPE TELESCOPES 


For a ray-trace analysis one must have the radii, eccentricities, and 
diameters of the mirror surfaces as well as the mirror separations. The 
eccentricities are often not readily available and are not easily measured. 
This appendix outlines the formalism used to describe the mirror surfaces foi 
the ray trace calculation and to show how the eccentricities are calculated. 
The approach follows the development of Jones (ref. 4) and can be applied to 
any Cassegrain- type telescope with reflectors made of conic sections. The 
classical Cassegrain and Dall-Kirkham geometries are treated explicitly. 

A conic section with a vertex at the origin and symmetric for x > O 1 is 
given by 

y 2 - 2rx + (1 - e 2 )x 2 = 0 (Al) 


where e >_ 0 is eccentricity, r = ed, and d is the distance from the focus 
to the directrix. The distance from the vertex to the nearest focus of the 
conic is 


x F = 


r 

1 + e 


(A2) 


The four standard cases are 


Figure 

Usual formula 

e 

r 

Circle 

Ellipse 

x 2 + y 2 = r 2 
x 2 /a 2 + y 2 /b 2 = 1 

0 

/l - b 2 /a 2 < 1 

r 

b 2 / a 

Parabola 

y 2 - 4ax = 0 

1 

2a 

Hyperbola 

x 2 /a 2 - y 2 /b 2 = 1 

A + b 2 /a 2 > 1 

b 2 /a 


For the ellipse and hyperbola, the distance between foci is 


2ae 


2 re 



(A3) 


From equation (Al) , x is an even function of y; with p = 1 - e 2 , one 
finds that 


1 A11 calculation in this appendix is done in a two-dimensional coordinate 

system with axes labeled x and y. 
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X 


■tH-sr 



8r 3 


+ sL 

16r 5 


r* + 


(A4) 


Consideration of equation (Al) for e = 0 clearly shows that r is the 
radius of the sphere that fits all the conic surfaces near the vertex - the 
so-called paraxial sphere. A conic section with its vertex at the origin may 
thus be characterized by parameters r, e, and maximum diameter D. 

In a classical Cassegrain, the primary surface is a paraboloid and the 
secondary is a hyperboloid. The secondary is chosen so that one of its foci 
coincides with the prime focus and the other is the Cassegrain focus. In this 
geometry, an incident beam of light parallel to the axis focuses to a point at 
the Cassegrain focus. The determination of the mirror parameters is straight- 
forward: the effective primary diameter Dp and focal length Fp are usu- 

ally ordained by practical considerations. Of course, ep = 1 and rp = 2Fp. 

The distance L from the prime focus to the Cassegrain focus is given by 
equation (A3) : 


L = 


2r s e s 
1 - e s 2 


(A5) 


The plate scale of the telescope is the ratio of the object angular size to 
the image physical size and is determined by the equivalent focal length Fp'; 
this is the focal length of a parabola of diameter Dp which has the same 
cone angle as the secondary: 


F 


P 



(A6) 


Here F s is the distance from the secondary to the Cassegrain focus, D s is 
the effective diameter of the secondary, and A is the amplification. 

Sketch (a) shows the various dimensions. (In the notation used in the text. 



Sketch (a) 
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The quantity D s satisfies the geometrical relation 



From equations (A6) and (A7) , L - F s = L/ (A + 1) so that 


(A7) 


D S 


L 

(A + 1) (Fp/Dp) 


(A8) 


Note that F p /D p is the f number of the primary and D s need not be the 
physical diameter of the secondary. (If the secondary is undersized to avoid 
filling the primary, then D s is the physical diameter of the secondary and 
Dp is the effective primary diameter.) Of course, equation (A8) is independ- 
ent of the eccentricities of the surfaces. From equation (A2) and sketch (a). 


L - F s 


1 + e< 


The combination of equations (A5) through (A7) and (A9) yields 


(A9) 


r c = 


eo = 


2 LA 

(AlOa) 

A 2 - 1 

A + 1 

A - 1 

(AlOb) 


for the parameters of the secondary in a true Cassegrain system. The result 
for r s is the same if one uses the spherical mirror formula, equation (2), 
which demonstrates the advantage of this approach, namely, that the funda- 
mental mirror parameters are the eccentricities of the conics and the radii of 
the corresponding paraxial spheres . 

Consider now the modification of the Cassegrain geometry to include other 
(conic) surfaces. Since equations (A8) and (AlOa) still hold, one needs to 
determine only e p and e s . According to equation (A4) , one can write for the 
primary surface 


y p 2 

X P 2% 


Zjl 


(2r p ) 


and, for the secondary surface. 


x 


s 




4A y s 4 
(A - 1)2 J (2r s ) 3 


(Alla) 


+ 


(Allb) 



Here P and S are parameters that describe the deformed surfaces (P * S * 0 
for the true Cassegrain geometry) and equation (AlOb) has been used to elimi- 
nate the eccentricity of the hyperboloidal secondary. The coordinates of each 
surface are measured from its intersection with the x axis. Deviations from 
the standard surfaces are 



(A12a) 


(A12b) 


for the primary and secondary, respectively. Spherical aberration can be 
corrected if the deformations of the primary and secondary surfaces are equal 
for corresponding zones, that is. 


with 


Ax p (y p ) = Ax s (y s ) 


(A13a) 


Zr-Sl 

y s D s 


(A13b) 


From equations (A12) and (A13), the condition for zero spherical aberration is 


- ■ *©■©■ 

Equation (A14) thus relates the deformation parameters for the surfaces of a 
modified Cassegrain telescope corrected for spherical aberration. 

For example, the Dall-Kirkham geometry uses a spherical secondary to 
facilitate construction and modification of the telescope. Comparing equa- 
tions (A4) and (Allb), one finds that 


S = 


1 + 4A - 
(A - l) 2 


(A15) 


since e s = 0. From equations (A14) and (A15) and the fact that P = 1 - e p 2 , 
it follows that 


■ [‘ - (^i)W(!)T /2< 1 


(A16) 
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so that the primary surface is ellipsoidal. If expression (AlOa) is 
substituted for r s and (A8) for D s into (A16) , one obtains, for the 
Dall-Kirkham geometry, a primary eccentricity of 


e P = 


(■ 


A 2 - 1 
A 3 



(A17) 
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APPENDIX B 


CALCULATION OF RESTRICTIONS ON USEFUL APERTURE 


Equation (5) , Dpu = FpD s /(Fp - z s ) , gives the underfilled primary 
diameter Dp u for an on-axis point source if one assumes that the secondary 
mirror is intentionally undersized. Here approximate expressions are given 
for increments to (l/2)Dp u due to other common effects. These expressions 
are useful in determining the secondary mirror diameter (D s ) when one designs 
the secondary and in determining the amount of the primary used for a given 
application once the telescope has been configured. The notation used here is 
the same as that in the text. 

A detector with a finite aperture extending a distance x d from the 
telescope axis in the image plane will detect light reflected from a point on 
the primary a distance (1/2) Dp U + Axx from the axis where 


Ax x « 


*d 2 s 
z s - z L 


(Bl) 


If the detector aperture is chosen to diffraction- limit the telescope at 
wavelengths greater than A d , then equation (9) applies with A d = A c and 
x d = (1/2 ) d . From equations (9) and (1), equation (Bl) becomes 


Axj = 


1.22A d z s 


(B2) 


To mount different instruments on the telescope, some adjustment of the 
focus must be provided. The primary aperture will increase as the image plane 
is moved closer to the primary for a secondary of fixed diameter. If Azj, is 
the distance the image plane is moved from the nominal position, the increase 
in the radius of the primary aperture is 


Ax 2 ~ 


D pu Az L 

(F p - z s )6 


(B3) 


where 6 is given by equation (4). 

Equations (Bl) to (B3) assume that the secondary is untilted. These 
effects are not included in the ray trace calculation. When the secondary is 
tilted to center the image of an object which is off-axis by an angle y, the 
extent of the primary used in the plane of the tilt is AX3 + (l/2)Dp u , where 

ZcFvj 

Ax 3 < \ y (B4) 

tp - z s 
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This effect is calculated more accurately by the ray trace program. The value 
obtained from equation (B4) is an overestimate. 

As an example, consider the Ax terms corresponding to the first case in 
table 1. If a 5-mm-diameter detector aperture on axis (corresponding to an 
8.85 arc-min full -width FOV) is assumed, then x^ = 2.5 mm and 


Axi = 1.6 mm 


(B5a) 


If = 150y, then 


Axi 1 = 0. 8 mm 


(B5b) 


Note that this corresponds to a diffraction diameter d = 2.5 mm. If the image 
plane is moved closer to the primary by Az^ = 25 mm, equation (B3) gives 


Ax 2 - 7.5 mm 


(B5c) 


For a beam throw of y = 5 arc-min. 


AX 3 5 3.3 mm 


(B5d) 


The true value of AX 3 obtained from the ray trace program (figs. 2 (b) and 
2(c)) is 1.9 mm. Since, for this example, Dp U = 289 mm and Dp = 305 mm, 

Axi, Axj ' , and Ax 3 are not sufficient to fill the primary. On the other hand, 
the value of Ax 2 would cause the primary to be filled and presumably be 
unacceptable for work beyond about 3y, due to the power reaching the detector 
from the telescope support structure. 
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APPENDIX C 


DETAILS OF RAY TRACE CALCULATION 


Rays Lj, L 2 , and L3 have direction cosines c^, c 2 j, and C3j, 
respectively, with j = x, y, and z. From appendix A, the equation of the 
primary is written 


$p = x 2 + y 2 + Cl - ep 2 )z 2 - 2rpZ = 0 


(Cl) 


Since secondary pivots through angle 0 about an axis parallel to the x 
axis through the pivot point at E, the equation of the secondary is 


where 


and 


$5 = x 2 + y' 2 + (1 - e s 2 )z' 2 - 2r s z’ = 0 ^ 

y' = (z - z 0 )sin 0 + (y - y o )cos 0 , 

z' = (z - Zq) cos 0 - (y - y Q )sin 9 , 

y Q = c sin 0 

z 0 = z s + c(l - cos 0) 


(C2) 


Here c 
which is 
the form 


is the distance from E to the center of the secondary surface, 
located at y Q , z Q . Equation (C2) can conveniently be rewritten in 


where 


$ s = x 2 + a! (y - y G ) 2 + a 2 (z - z 0 ) 2 + 2a 3 (y - y 0 ) (z - z Q ) 
+2a 4 (z - Zq) + 2a 5 (y - y Q ) =0 


a i = 1 - e s 2 

sin 2 

a 2 = l - e s 2 

cos 2 

a 3 ~ e s 2 sin 20 

a 4 = -r s cos 

0 

ag - r s sin 0 



(C3) 


The unit outward normal to one of these surfaces is 

~ V$i 


Wi\ 


(C4) 
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where i = p or s. Given Li , one has 


Lj • h p = -L 2 • n = cos a > 0 \ 

/ } (C5) 

|Li X npj = |L 2 x n p | = sin a >_ 0 J 

The equation of L 2 is then obtained from equation (B5) : 

£.2 = Li - 2rip cos a (C6) 


If 



then points lying on L 2 satisfy 


(C7) 


x - x A — m 2x (z z A ) 
y - y A = m 2y (z - z A ) 


CCS) 


where A denotes the intersection of L x with the primary. 

The simultaneous solution of equations (C3) and (C8) gives the intersec- 
tion of L 2 with the secondary at point B. One finds a quadratic equation 
for u = zg - z A , with the result 


u 


B + (B 2 - AC) 1 / 2 
A 


CC9) 


where 

2 2 

A = m 2x + aim 2y + a 2 + 2a 3 m 2y 

B = m 2x x A + aim 2y w y + a 2 w 2 + a 3 (w y + m^vi z ) + a 4 + a 5 m 2y 
C = x A 2 + aiw y 2 + a 2 w z 2 + 2a 3 w y w z + 2a 4 w z + 2a 5 w y 
w j " ^A " jo » j = x, y, z 


(CIO) 


The coordinates of point B are then 

x B = x A + m 2x u 
Xb = X A + m 2y u 


CC11) 
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z B = z A + u 



The same procedure which led to equation (C 6 ) gives 


^ A A 

L 3 = L 2 + 2 n s cos g 


where equations (C3) , (C4) , and (Cll) are used to obtain h s and 
cos g = -£2 * n s . As in equation (C7) , one defines 


m 3 



If the image plane is at z = zl> then the intersection with £3 is 
given by 


X D “ X B + ( Z L " ZB^x' 
Xd = ye + ( Z L - z B) m 3 y ( 


which is the desired result. 


(C12) 


(C13) 


(C14) 
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